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1. Discovery of superfluid cooling of the Cassiopeia A neutron star 

Neutron stars are created in the collapse and supernova explosion of massive stars, and they 
begin their lives very hot (with kT > \Q MeV) but cool rapidly through the emission of neutrinos. 
This neutrino emission depends on uncertain physics at the supra-nuclear densities (> 0.08 fm^^) 
of the neutron star core [jlj ^ ||, Current theories indicate that the stellar core may contain ex- 
otica, such as hyperons and deconfined quarks, and matter may be in a superfluid/superconducting 
state [^, ^ ^. By observing the cooling of neutron stars and comparing their temperatures to 
theoretical models, we can constrain the nuclear physics properties that govern the stellar interior. 

The compact object at the center of remnant of the Cassiopeia A supernova was discovered 
in Chandra X-ray Observatory first-light observations and subsequently identified as a neutron 
star [^. The supernova explosion is estimated to have occurred in the year 1681 ± 19 [|l^]; this 
makes the Cassiopeia A neutron star the youngest-known neutron star at an age of 330 yr. A 
steady temperature decline of four percent was found using Chandra observations taken during the 



last 10 years [11]. If the rapid decline is due to passive neutrino cooling, then this is the first direct 



evidence for superfluidity and superconductivity in the core of a neutron star [12, [131 ]. 

The left panel of Fig. |l] (in particular, see inset) shows Chandra temperature measurements of 
the Cassiopeia A neutron star from 1999 to 2010 [ pT| |l3]|. Figure |T]also shows surface temperatures 
for three theoretical models of neutron star cooling: "N — normal matter" corresponds to neutron 
star matter that does not contain any sort of superfluid, "pSF — proton superfluid" is for superfluid 
protons in the core, and "npSF — neutron/proton superfluid" is for superfluid neutrons and protons 
in the core. Note the difference between the cooling behavior of models with normal matter (N) and 
matter containing supeiiluids (pSF or npSF) after 40 yr. In the latter models, a proton supercon- 
ductor foiTiis soon after neutron star formation, and this suppresses neutrino emission, so that the 
cooling rate is weaker than for normal matter. This enables the star to stay relatively wami, leading 
to a rapid temperature drop once neutrons become superfluid [|T^, 15]. The model with superfluid 



neutrons and protons (npSF) fits the data at an age of a few hundred years. The four circles trace 
the cooling curve predicted by this model from about 10 years after the supernova explosion (SN 
in ^ 1680) to about the time when neutrons become superfluid in the core: (I) At early ages, the 
neutron star core cools so rapidly by neutrino emission that the crust does not have time to react. 
Thus the crust is hotter than the core in 1690 (age k, 10 yr; protons are superconducting by this 
time), and the surface temperature declines very slowly. (2) The surface temperature eventually 
reacts to the "cooling wave" that sweeps through the crust and starts to drop off more quickly. 
After 1760, the temperature becomes almost constant throughout the star. (3) Then in « 1900, the 
interior temperature drops below the critical value for a neutron superfluid to form and enhanced 
neutrino emission occurs in the core, as neutron Cooper pairs form. Energy is lost as the neutrinos 
are emitted, causing the core to cool off and another cooling wave to travel outwards. As neutrons 
in large regions of the core become superfluid, the surface temperature drops off quickly, beginning 
in « 1930 (i.e., start of "Great Depression") and continuing through the present date. See Fig. 2 of 



[]16|], which shows evolution of interior temperature T{p) and transition to neutron supeiiluidity. 

Monitoring of the temperature decline will allow improved constraints on the (1) critical tem- 
perature for neutron triplet pairing Tnt (maximum pairing gap energy; see Fig. |I| and Fig. 1 of [13]), 
(2) suppression due to collective effects of the axial vector current for pair formation (see Fig. 3 of 
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Figure 1: Left: Theoretical models of neutron star cooling with superfluid neutrons and protons (npSF — 
solid), normal neutrons and superfluid protons (pSF — long-dashed), and normal neutrons and protons (N 
— short-dashed). Circles indicate temperature for the npSF model at particular times/years. Crosses are 
Chandra X-ray Observatory measurements of the Cassiopeia A neutron star Right: Superfluid pairing gap 
energies as a function of Fermi wavenumber for neutrons and protons k^^. The maximum neutron triplet 
gap is taken to be either the shallow model or the deep model, and each triplet gap leads to cooling that can 
fit the Cassiopeia A data (see [ [T2[ | and [[T3[], respectively). 



[13]), and (3) neutron star mass and nuclear equation of state (see Fig. 1 of [|Tj] and Fig. 4 of [12]). 

Guided by the discovery of a superfluid and superconductor in the Cassiopeia A neutron star, 
we examine three examples where measurement of the pairing gap energies has possible effects or 
where further constraints may be obtained. 



2. Nuclear X-ray bursts and neutrino cooling by neutron superfluid 

In contrast to the neutron star in Cassiopeia A, many old neutron stars are found in binary 
systems. These binaries can be seen in X-rays, which are produced when material from the com- 
panion star accretes onto the neutron star. If the companion has a low mass, the systems are known 
as low-mass X-ray binaries (LMXBs). Many LMXBs undergo bright X-ray bursts due to unstable 
thermonuclear burning of hydrogen and/or helium in the surface layers of the neutron star [17, p^ ]. 



Bursts are sometimes observed to recur in individual sources, and recurrence times between multi- 



ple bursts span a wide range, from minutes to days QlSj |19| ]. However, recurrence times < 1 hr are 



too short for the neutron star to accrete enough fuel for subsequent bursts [|T7[ 



In [|2ip, we revisit the method used to infer core temperatures Tc of neutron stars in LMXBs. 
Compression by accreted matter induces nuclear reactions in the deep crust, which release 1.5 
or 1.9 MeV nucleon^^ [22], and this heats the core directly by a luminosity Lheat ~ 0.0078 L, 



aco 



where Lace is the time-averaged X-ray luminosity of the LMXB [ ^3} p4| ]. Figure ^ shows the 
measured heating rate Lheat > as well as the theoretical neutrino luminosity Ly, which depends on 
the neutron tiiplet pairing gap energy (or critical temperature T^x). The intersection of the curves 
Lheat and Ly yields the neutron star core temperature. We see that the core temperature in relatively 
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high luminosity LMXBs may not be uniquely determined. If Tnt^max ^ 8 x 10^ K, there can be 
two thermally stable values' of the core temperature associated with a single observed accretion 
luminosity. For example, there is a factor of < 3 difference in the inferred if Lace ~ (0.2 — 
9) X 10^'' ergs s"' and rnt,max = 4.3 x 10^ K. The luminosities of all LMXBs that show short burst 
recurrence time lie within this range [19|. To highlight this point, we place the LMXBs with "short 
time" on the high-temperature branch and LMXBs with "long time" on the low-temperature branch. 
Thus the sources with short time and higher temperatures have normal neutrons in the stellar core, 
while sources with long time and lower temperatures have superfluid neutrons. 

If short burst recurrence time LMXBs do indeed possess hotter core temperatures, then mea- 
surements of the minimum and maximum accretion luminosities of bursts from short time LMXBs 
and long time LMXBs, respectively, can be used to constrain the neutron superfluid critical temper- 
ature Tni{p). This is illustrated in Fig. ^ where it is clear that the accretion luminosities for LMXBs 
can constrain Tnt.max and Pnt,peak> while the width of Tai{p) is not as important in determining the 
qualitative behavior of Ly. 




Figure 2: Left: Simple models of the neutron triplet critical temperature Tlit(p). Right: Neutrino luminosity 
as a function of neutron star core temperature, where the different curves are Ly calculated using the models 
of Tni{p) shown in the left panel. Upper (lower) horizontal dot-dashed line is the highest (lowest) observed 
Lheat from among all long (short) recurrence time bursts. Squares and diamonds are where Lheat = (with 
4.3 X 10^ K,9.4 X lO''* g cm^^) for each LMXB with long and short recurrence times, respectively. 



3. Neutron star spin and damping of r-mode oscillations 

One of the main mechanisms that is expected to affect the spin evolution of an accreting 
neutron star in a LMXB is the instability associated with r-modes, which are a class of oscillations 
in a star whose restoring force is the Coriolis force. The emission of gravitational waves can 



excite r-modes in the stellar core and cause the oscillations to grow [25]. The r-mode instability 



is interesting for many reasons, mainly because the associated gravitational wave signal may be 



'There can be three values of Tq that intersect each horizontal Lheat- But the intermediate temperature is thermally 
unstable since a temperature decrease leads to an increase in neutrino luminosity which causes even more rapid cooling. 
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detectable, but also because its understanding requires knowledge from a wide range of nuclear 
physics. A primary agent that enters the r-mode discussion is damping mechanisms related to 



shear and bulk viscosities and exotica like hyperons, quarks, and superfluid vortices [ |26| ] (see also 
C.J. Horowitz, this volume). The instability depends primarily on the neutron star spin rate Vg and 
core temperature Tc. This leads to an instability "window," determined by a critical curve in the 
V^-Tc plane, inside which the instability is active. What has not been appreciated is that this leaves 
the majority of the observed systems significantly inside the instability window: rapidly rotating 
neutron stars (i.e., pulsai^s) should not possess spin rates at their observed levels [^, 

One solution to this dilemma is to change the window so that r-mode growth is stabilized at 
relatively high spin rates. However in order to do this, a revision of our understanding of damping 
mechanisms is required. In [^], we explore the possibilities. For example, there may be resonances 
between the r-mode and torsional oscillations of the elastic crust [p9|]. Such resonances could have 
a sizeable effect on the instability window. Figure ^ shows an example; the illustrated instability 
window has a broad resonance at 600 Hz, which is the typical frequency of the first overtone of 
pure crustal modes. 

Another possibility is an instability spin frequency that increases with temperature in the range 



of interest here [ pO[ , [31| |. If this is the case, then neutron stars may evolve to a quasi-equilibrium 
where the r-mode instability is balanced (on average) by accretion and r-mode heating is balanced 
by cooling. This solution is interesting because it predicts persistent (low-level) gravitational ra- 
diation. Figure ^ shows a model using hyperon bulk viscosity suppressed by superfluidity. This 
explanation has a major problem though. It must be able to explain how observed pulsars with 
millisecond spin rate emerge from accreting systems. Once the accretion phase ends, the neutron 
star will cool, enter the instability window, and spin down to ~ 300 Hz. In other words, it would 



be very difficult to explain the formation of a pulsar spinning at 716 Hz [32]. 

An intriguing possibility involves mutual friction due to vortices in a rotating superfluid. The 
standai^d mechanism (electrons scattered off of magnetized vortices) is too weak to affect the in- 



stability window [p3|]. However, if we increase (arbitrarily) the strength of this mechanism by a 
factor ~ 25, then mutual friction dominates the damping, as shown in Fig. |3[ Moreover this would 
set a spin threshold for instability similar to the highest observed Vg and would allow systems to 
remain rapidly rotating after accretion shuts off. Enhanced friction may result from the interaction 
between vortices and proton fluxtubes in the outer core. 

4. Pulsar glitches and neutron superfluidity in the crust 

Mature neutron stars tend to have extremely stable spin rates, with some pulsars possessing 
a timing stability that rivals the best terrestrial atomic clocks. However, young neutron stars may 
behave in a less ordered fashion. In particular, many young pulsars exhibit regular glitches, where 



the obsei-ved spin rate suddenly increases [36|. The consensus view is that these events are due to 



a superfluid component in the stellar interior [37]. Anderson & Itoh ]pq] envisaged a glitch as a 



tug-of-war between the tendency of the neutron superfluid to match the spindown rate of the rest of 
the star by expelling vortices and the impediment experienced by moving vortices that are pinned to 
crust nuclei. Strong vortex pinning prevents the superfluid from spinning down and creates a spin 
lag with respect to the rest of the star (which slows by magnetic dipole radiation). This situation 
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Figure 3: Three scenarios that could explain r-mode stability in the observed LMXBs (squares and dia- 
monds; see Fig. 2). R-mode growth is stable below (i.e., at lower Vs) the various curves, while the dashed 
lines at 900 Hz indicate the break-up limit. Left: Crust mode resonance at 600 Hz. Middle: Superfluid 
hyperons (based on [|4| with % = 01). Right: Strong vortex mutual friction (based on the strong/weak 
superfluidity models from [B5| with ^ w 0.01). 



cannot persist forever. The increasing spin lag leads to a build-up in the Magnus force exerted on 
the vortices. Above a threshold, pinning can no longer be sustained, vortices break free, and excess 
angular momentum is transferred to the crust. This leads to the observed spin-up, i.e., glitch. 

A previous analysis by [39] suggests that glitches involve a superfluid reservoir with moment 
of inertia /n// ~ 1%, where / and /„ are the moments of inertia of the entire star and the neutron 
superfluid component, respectively. The similarity of the infeixed /„ to the theoretically estimated 
moment of inertia of the crust (which is dominated by free neutrons in the inner crust) for realistic 



nuclear equations of state [ |40| ] supports the idea that glitches involve only the crust region. In [pHQ, 
we show that this logic breaks down when one accounts for non-dissipative entrainment coupling 
between the neutron superfluid and the crust lattice, an effect which can be expressed in teniis of 
an effective neutron mass m*. Recent work indicates that this effective mass may be significantly 
larger than the bare neutron mass «!„ [42] (see Fig. This implies a decreased superfluid mobil- 
ity with respect to the lattice and the need for a lai^ger angular momentum reservoir for glitches. 
Combining the latest glitch data [36] with a general relativistic multifluid model that includes en- 
trainment, we find that the requisite superfluid moment of inertia is above the capacity of the crust 
superfluid [41] (see Fig. ^. Some solutions are briefly discussed below (see also [pd]]). 

One possible explanation could be that the superfluid in the core is involved in the glitch (see, 
e.g., [44]), and the combined superfluid moment of inertia reservoir is just large enough to explain 
the observations. If this fine-tuning resolves the problem, then a more detailed calculation would 
constrain the singlet pairing gap for neutrons. This would be an interesting complement to the 
constraints on core superfluids (singlet protons and triplet neutrons) discussed in Sections |]-|3[ 

Another solution could be the result of superfluid behavior at the cmst-core transition. Unless 
the superfluid is confined to the crust, one would have to explain why the crust component decou- 
ples from the core during the glitch event. This would be particularly vexing if the singlet pairing 
gap is such that the neutron superfluid reaches far into the core. A central issue concerns the na- 
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Figure 4: Left: Neutron effective mass as a function of baryon number density. Triangles are from [M 



while the curve is a fit from |43|. Right: Moment of inertia ratio /„// as a function of mass for neutron 
star models built from the APR I (solid) and SLy (dotted) nuclear equations of state. If glitches in the Vela 
pulsar are to be explained solely by a crust superfluid, then the moment of inertia ratio must satisfy /„// > 
0.016 X (< m* > /nin) k, 0.07, where the average effective mass is calculated from that shown in the left 
panel; also shown is the constraint when entrainment is not taken into account, i.e., when < m* > /m^ = 1. 



ture of superfluid vortices extending across this interface. The standard picture is that vortices are 



magnetized in the core [|45|], due to entrainment and the presence of superconducting protons, but 
not in the crust. This suggests a more compUcated transition behavior than is usually assumed. 

Full details of the work presented in Sections |l]-|| can be found in [|l|] (see also [|l|]), [||], 
[|7|] (see also and @ (see also [^), respectively. 



WCGH is indebted to Daniel Patnaude, Peter Shtemin, and Dmitry Yakovlev for assistance. 
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